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Abstract 



We have derived the local electric field inside graded metal-dielectric com- 
posite films with weak nonlinearity analytically, which further yields the ef- 
fective linear dielectric constant and third-order nonlinear susceptibility of 
the graded structures. As a result, the composition-dependent gradation can 
produce a broad resonant plasmon band in the optical region, resulting in 
a large enhancement of the optical nonlinearity and hence a large figure of 
merit. 
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In contrast to bulk materials, the corresponding thin films often possess different optical 
properties [1,2]. It is also known that graded materials [3,4] have quite different physical 
properties from the homogeneous materials [5]. To one's interest, one found that the graded 
thin films may have better dielectric properties than a single- layer film [6] . However, the tra- 
ditional theories [7,8] fail to deal with the composites of graded particles. For this purpose, 
we put forth a first- principles approach [9,10] and a differential effective dipole approxima- 
tion [11]. The situation becomes more complicated by the presence of nonlinearity which 
exists in real composites [12-22] . Moreover, there is a great need for nonlinear optical mate- 
rials with large nonlinear susceptibility or optimal figure of merit (FOM). Therefore, much 
work has been done on how to gain a large nonlinearity enhancement or optimal FOM of 
bulk composites by the surface-plasmon resonance in metal-dielectric composites [19], as well 
as by taking into account the structural information [23,24]. Recently, a large nonlinearity 
enhancement was found for a sub-wavelength multilayer of titanium dioxide and conjugated 
polymer [25]. 

However, the surface plasmon resonant nonlinearity enhancement often occurs concomi- 
tantly with a strong absorption, and unfortunately this behavior renders the FOM of the 
resonant enhancement peak to be too small to be useful. To circumvent this problem, we 
shall consider a kind of graded metal-dielectric composite film, in which a dielectric compo- 
nent is introduced as spherical particles embedded in the metallic component. 

Let us start by considering a graded film with width L. Here the gradation of interest is 
in the direction perpendicular to the film. In this connection, the local constitutive relation 
between the displacement D and electric field E is given by 

ID{z,uj)^e{z,uj)E{z,uj)+x{z,uj)\E{z,uj)\^E{z,uj), (1) 

where e{z,u!) and xi^,^) are respectively the linear dielectric constant and third-order non- 
linear susceptibility of a certain layer inside the graded film. It is worth noting that both 
e{z, uj) and x(2;, ui) are gradation profiles as a function of position z. In view of the metal- 
dielectric composite layer, e{z, lu) is given by the well-known Maxwell-Garnett approximation 
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62 - ei{Lj) 
62 + 2ei(a;) ' 



(2) 



e{z,uj) + 2ei{uj) 



where 62 stands for the dielectric constant of the dielectric, and p{z) the layer dielectric 
profile. It is worth noting that Eq. (2) gives an approximate e{z, uj) for the present geometry. 
Here the dielectric function of the metal €1(00) is given by the Drude expression 



where Up denotes the bulk plasmon frequency, and 7 the collision frequency. Regarding 
Eq. (2), we should remark more. In fact, it is not possible to calculate e{z,u) exactly in 
terms of the layer dielectric profile p{z). However, to obtain an estimate of e{z,u!), we can 
take a small volume element inside the layer, at a position z. Further, this small volume 
element can be seen as a composite where the dielectric particles are randomly embedded 
in the metallic component. Accordingly, the volume fraction of the dielectric particles is 
p{z). In this regard, the Maxwell-Garnett approximation [namely, Eq. (2)] holds well for 
computing e{z,uj). 

Let us further assume that the weak nonlinearity condition is satisfied in the present work. 
That is, the contribution of the second term (nonlinear part x(2;, a;)|E(2;, a;)p) in the right- 
hand side of Eq. (1) is much less than that of the first term (hnear part e{z,u)) [12]. Next, 
we focus on the quasi-static approximation, under which the whole graded structure can be 
regarded as an effective homogeneous one with effective (overall) linear dielectric constant 
e{cu) and effective (overall) third-order nonlinear susceptibility x{^)- the mathematical 
expression, the definition of e(a;) and xi'^) is given by [12] 



where (• • •) stands for the spatial average of • • •, and Eq = EqCz the applied field along 
axis. 

Owing to the simple graded structure, we can use the equivalent capacitance of series 
combination to calculate the linear response (i.e., optical absorption for the metallic film). 



eiico) = 1 




(3) 



(D) = e(a;)Eo + x(^)|EorEo, 



(4) 
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For calculating the nonlinear optical response, we first calculate local electric field E{z, u) 
by the identity 

e{z,u!)E{z,u!) = e{u!)EQ (6) 

due to the virtue of the continuity of the electric displacement. In view of the existence 
of nonlinearity inside the graded film, the effective nonlinear response xi^) can be given 
by [12] 

xHEo^ = {x{z,u;)\E,^{z)\'E,^{zf), (7) 

where E^n denotes the linear local electric field. We can take one step forward to express 
the effective nonlinear response as an integral over the film. 



1 



LJo 



e{z,u;) 



In fact, the assumption of a 2;-dependent Maxwell-Garnett type dielectric profile is equivalent 
to assume that the local field depends on z only, see Eq. (2). However, the real x(a;) should 
involve an integral over x, y, and z of the local xi^j Ui multiplied by terms involving 
e{x,y, ZjLu). Thus, Eq. (8) ofi^ers an approximate x(to'), as expected. 

In what follows, we shall do some numerical calculations. First, set xi^^^) be a 
constant Xi- By doing so, the two components are assumed to have the same real and 
positive frequency- independent Xi) so that we could emphasize the enhancement of the 
optical nonlinearity. Regarding the layer dielectric profile, let us take a power former p{z) — 
az^. Without loss of generality, the layer width L is taken to be 1. 

In Fig. 1, we display (a) the optical absorption ~ Im[e(c(j)], (b) the modulus of 
the effective third-order optical nonlinearity enhancement |x('^)|/Xi) ^^'^ is) the FOM 
|x(<^)|/{XiIn^[^('^)]} a function of the incident angular frequency uj. Here Im[- • •] means 
the imaginary part of • • • . To one's interest, when the layer dielectric profile p{z) is taken 



into account, a broad resonant plasmon band is observed always. In other words, the broad 
band is caused to appear by the effect of the positional dependence of the dielectric or metal. 
Also, we find that increasing a causes the resonant band not only to be enhanced, but also 
red-shifted (namely, located at a lower frequency region). In this work a should satisfy 
< a < 1, in order to ensure a p{z) which is always smaller than unity. In a word, although 
the enhancement of the effective third-order optical nonlinearity is often accompanied with 
the appearance of the optical absorption, the FOM is still possible to be very attractive 
due to the presence of the positional dependence of the dielectric or metallic components. 
Moreover, it is worth noting that a prominent surface plasmon resonant peak appears at 
somewhat higher frequencies in addition to the surface plasmon band. As a increases, this 
peak is blue-shifted (i.e., locates at a higher frequency region) accordingly. 

Similarly, Figure 2 displays the influence of m. It is apparent to see that the broad 
resonant plasmon band can be enhanced significantly by adjusting m. However, no distinct 
red-shift occurs for the plasmon band as m varies. In contrast, we notice that increasing m 
can make the surface plasmon resonant peak red-shifted. 

As a matter of fact, the present results do not depend crucially on the particular form of 
the layer dielectric profile p{z) . The only requirement is that we must have a composition- 
dependent layer to yield a broad plasmon band for the graded film. It should be remarked 
that the optical response of the graded structure depends on polarization of the incident light, 
because the incident optical field can always be resolved into two polarizations. However, a 
large nonlinearity enhancement occurs only when the electric field is parallel to the direction 
of the gradient [25] , and the other polarization does not produce nonlinearity enhancement 
at all [25] . Fortunately, the nonlinear susceptibilities of both the parallel and perpendicular 
polarizations are related to the nonlinear phase shift, which can be measured by using the z- 
scan method [25] . It is of interest to extend our consideration to composites in which graded 
spherical particles are embedded in a host medium [26] to account for mutual interactions 
among graded particles [24]. 

To calculate e{z,u!), we used Eq. (2) in which metal (or dielectric) serves as a host (or 
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inclusion). Inversely, we could see metal (or dielectric) as an inclusion (or host), and use 
the same form as Eq. (2) by exchanging ei{uj) and €2, and p{z) and 1 — p{z). In case of 
ei(a;) > 62, the former always gives the upper bound, while the latter the lower bound, and 
vice versa. The exact result must lie between the two bounds [27] . 

To sum up, we have studied a graded metal-dielectric composite film, and found that 
the composition-dependent gradation in the film yields a broad resonant plasmon band in 
the optical region. Thus, it is possible to gain a large nonlinearity enhancement as well as 
optimal FOM for graded metal-dielectric composite films. 
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FIGURES 

FIG. 1. (a) Linear optical absorption Im[e(a;)], (b) enhancement of the third-order optical 
nonlinearity |x('^)l/xi) (c) FOM (figure of merit) = |x(<^)|/{xil™[^('^)]} versus the normalized 
incident angular frequency co/up for layer dielectric profile p{z) = az^, for different a at m = 1.0. 
Solid line: a = 0.2; Dashed line: a = 0.4; Long-dashed line: a = 0.6; Dot-dashed line a = 0.8. 
Parameters: ^/ujp = 0.01 and 62 = (3/2)^. 

FIG. 2. Same as Fig.l, but for different m at a = 0.8. Solid line: m = 0.2; Dotted line: m = 0.6; 
Dashed line: m = 1.0; Long-dashed line: m = 1.4; Dot-dashed line m = 1.8. 
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